function col_of_plotIntensity(T, r, z, time_step)
    % 获取对应时间步的温度场
    T_at_time = T(:,:,time_step);
    T_at_time1 = T_at_time / max(T_at_time(:)) * 100;  % 归一化温度场

    % 获取第500行的数据（r方向上的强度）
    intensity_data = T_at_time1(:, 500);

    % 绘制强度图
    figure;
    plot(r, intensity_data, 'LineWidth', 2);
    xlabel('r (mm)');
    ylabel('Temperature (%)');
    title(['Temperature Intensity at Row 500 (t = ', num2str(time_step), ' s)']);
    grid on;
end
